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g ■ Abstract 

O. 

^P , The Wheeler-DeWitt equation is solved for the Bergmann- Wagoner scalar- 

tensor gravitational theory in the case of Friedmann-Robertson- Walker cos- 
mological model. We present solutions for several cosmological functions: 
i) X{(p) = 0, a) X{(f)) = 3Ao(/> and Hi) a more complex A(0), that depends 



o 



bJQ. on the choice of the coupling function, considering closed, flat and hyperbolic 

Friedmann universes (/c = 1,0,-1). In the first two cases we show parti- 
cular quantum wormhole solutions. Also, classical solutions are considered 



>< 



$H ' for some scalar-tensor theories, and we study the third quantization of some 

minisuperspace models. 
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I. INTRODUCTION 

In this work we will consider the quantum cosmology of the Bergmann- Wagoner (BW) 
theory of gravitation |]I| in the homogeneous and isotropic case. This theory is the most 
general scalar-tensor theory of gravitation. The action for this theory contains two arbitrary 
functions of the scalar field, uj{(j)) and A(0). Every specific choice of these functions define 
a particular scalar-tensor theory. The function uj{(f)) is a coupling function. A dual role 
is played by the function A(0), it acts as a dynamical cosmological function that can give 
cosmological models with small cosmological constant at present times consistent with the 
observations but large at earlier times to produce inflation, and is also a potential function 
in the equation for the scalar field. Among the particular cases of BW theory we have: 
(1) Brans-Dicke theory (a;(0) = cuq = const ., \{(f)) = 0)0, (2) low energy limit of string 
theory (c<;(0) = -1), Barker's constant G theory (u;(0) = (4 - 3(j))/{2(f) - 2), A(0) = 0)0, 
Kaluza-Klein theory {u = —{n — l)/n, where the space-time dimension is 4 + n) |Q. In 
the past several cosmological consequences of the classical BW theory have been studied: 
exact solutions, qualitative analysis of the equations of motion and the implementation of 
inflation ("no-hair" theorem). 

Very little attention has been paid to the quantum cosmology of this general theory, 
however several particular case has been studied 0. Therefore it seems of interest to con- 
sider the quantum cosmology obtained from the BW theory. In what follows we consider the 
quantum cosmology of BW theory, that is, we keep as far as possible co'(0) and A(0) arbi- 
trary in the homogeneous and isotropic line element, we obtain the Wheeler-DeWitt(WDW) 
equation p for this model. When we try to solve exactly the WDW equation by separation 
of variables, it is found that we can do so for arbitrary uj{(f)), but only for particular forms 
of A(0). We solve the WDW equation for several potentials, finding general solutions to 
the separated equations and from them we can use superposition to obtain wavefunctions 
that satisfy particular boundary conditions. In particular we consider the wormhole wave- 
functions that satisfy the HP regularity conditions [0 : (i) the wavefunction is exponentially 
damped for large spatial geometry, and (ii) the wavefunction should be well behaved (reg- 
ular) at the origin. This is because the wavefunction should represent Euclidean space for 
large radius and there should be no singularities when the spatial geometry degenerates. At 
the classical level, BW theory is known to have solutions representing wormholes [§ . Then, 
we explore the possibility to get wormholes in the quantum case. 

From the field equations one can also consider the classical problem; for those cases in 
which classical solutions were not known, they are given. Explicitely for Brans-Dicke theory. 
Barker's theory and three families of other theories that have been studied previously P,p!0|. 
In one particular case we obtained classical non-singular solutions. 

On the other hand, because the fact that the WDW equation is a second-order hyperbolic 
differential equation, there is no conserved positive-definite probability density as in the case 
of the Klein-Gordon equation. Then, there is a problem in the interpretation that | V' P 
is a probability, where ^ is a solution to the WDW equation. An alternative approach is 



the third quantization procedure |TT|, in analogy with the second quantization of the Klein- 
Gordon equation. Thus, it can be construct a consistent probability measure in quantum 
gravity by promoting ip to a. quantum field operator that acts on a Hilbert space of states. 
Then, third-quantized universe theory describes a system of many universes. For one of the 
potentials and the closed universe we carry out the third quantization of the model and 
calculated the number of universes produced from nothing. 

Another interesting aspect on the sense of third quantization is the study of quantum 
effects in the universe evolution, for a specific potential we calculated the uncertainty rela- 
tions and found that the quantum fiuctuations are large for small values of the expansion 
factor of the universe and are bounded for large values of it. 

The organization of the paper is as follows: in the Sections II-IV we give a short review 
of the minisuperspace models that lead to the class of WDW equations considered here, and 
respective quantum solutions are given. In Sec. V we show a collection of classical solutions 
for some scalar-tensor theories. In Sec. VI the quantum wormhole configurations are studied. 
Sec. VII is devoted to study the third quantization of particular minisuperspace model for 
N = l/a</)2. In Sec. VIII we calculate the fiuctuations of a third quantized minisuperspace 
model for A^ = 1. Finally, we summarize and conclude in Sec. IX. 



II. THE BERGMANN-WAGONER MINISUPERSPACE MODEL 

In order to construct the minisuperspace model, our starting point is the action for 
Bergmann- Wagoner scalar-tensor theory: 



•^=7^ / V-^ 



0i?{4)_i^^M^0^^0_^ + 20A(0) 



t/lr + 



dM 



Vh(j)hijK'^ d 



X, 



(2.1) 



where g = det{g^^^), R^^^ is the scalar curvature of the Friedmann- Robertson- Walker(FRW) 
theory, (j){t) is the conventional real scalar gravitational field, Ip is the Planck length and 
X{4>) is the cosmological term. The second integral is a surface term involving the induced 
metric hij and second fundamental form K^^ on the boundary, needed to cancel the second 
derivatives in R^^^ when the action is varied with the metric and scalar field, but not their 
normal derivatives, fixed on the boundary. We want to study an homogeneous and isotropic 
cosmological model, consequently we used the FRW metric line element in spherical polar 
coordinates (t, r, 9, $), given by 



ds' 



-N\t)dt^ + a^{t) 



dr"^ 



kr' 



+ r2(rf^2^sin2 0rf$^ 



(2.2) 



where A^ is the lapse function, a is the scale factor of the universe and the Ricci scalar is 



d{4) ^ _^ _ a_^ 

0? ma^ 
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(2.3) 



substituting Eq. ( |2.3| ) into Eq. ( |2.1|) and integrating with respect to space coordinates, after 
simplifications we have 
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dt. 



(2.4) 



where dot denotes the time derivative with respect to the time t, now we introduce the 
following new variables 



X = a(p2, y 



^^^ -f, A{y) = ^, dT = (j)-^dt, (2.5) 
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then the BW action simplifies to 



S 



BW 



—x'"^ - —y'^ - Nkx + Nx^A{y) 



dr. 



(2.6) 



prime denotes the time derivative with respect to the new time r. The lapse function A^ is a 
Lagrangian multiplier whose variation leads to a constraint equation. The Hamiltonian can 
be constructed according to the standard canonical rule, and we expressed it as iiT = NH. 
Since A^ is a Lagrange multiplier, we have the constraint ?i ~ 0, and follow the Dirac 
quantization procedure; the wavefunction of the universe is annihilated by the Hamiltonian 
operator HiIj = 0, this equation is know as WDW equation. We can see that this wave 
equation is independent of the lapse, however in the following sections we study particular 
gauges A^ = 1 and A^ = 1/x, that results in separable WDW equations. 

III. GAUGE N=l 



With the choice A^ = 1, the action ( p^.6| ) becomes 

xx'"^ — x^y''^ — kx + x'^A{y) 



Sbw — - 



dr. 



/ Ldr, 



-x^y'. 



hence the canonical conjugate momenta corresponding to x and y are given by 

_ aL _ , _dL_ 

ox' oy' 

The Hamiltonian H of the system is 

1 r 



n 



X TT^- X -^TT^ + kx - x'^A{y) 



(3.1) 



(3.2) 



(3.3) 



now, the canonical momenta in Eq. ( p.3| ) are converted into operators in the standard way, 
vr^ — > —x~^^ (•^^^j ^^"^ '^y ~^ ~W^i ^^^ ambiguity of factor ordering is encoded in the p 
parameter, then the WDW equation is 



(9^ (9 (9^ 

2 74 6a/\ 



7/'(x,2/) = 0. 



(3.4) 



dx"^ dx dy'^ 

It is easy to see that this wave equation can be solved by means of separation variables 
if A(|/) = const., thus in the next subsections we present three simple solvable cases for 
different cosmological terms. 



A. Case A(y) = 

First, let us consider the WDW equation without a cosmological term, then we have 
o 92 d d^ 



X 



px- 



kx 



i^{x,y) = 0, 



(3.5) 



dx"^ dx dy"^ 

in order to solve this wave equation we make ipix, y) = e*^^X(x), thus we obtain the following 
solution given in terms of Bessel functions. In the next, we show the corresponding universe 
wavefunction for the three possible values of spatial curvature constant k = —1, 0, 1. 

i) k=l 

In a closed FRW model we have that the universe wavefunction is 



1ps{x,y) = X 2''y2 



Cil„ 



X 

T 



CoK 



X 



2-f^m 



Cse'^y + de- 



isy 



(3.6) 



a) k = —1 

In the hyperbolic FRW universe, the solution is expressed in terms of modified first kind 
Bessel functions 

^s{x,y) = x"-^y"^ [Ci/^ {"C) + C2l-m dW fee^'^ + C^e 



■isy 



2 ) V 2 , 

where s is a separation constant, and the Bessel function order is given by m 



(3.7) 



S2. 



in) k = 

For a flat FRW universe the wavefunctions ip{x,y) = e'^^X(x) are 

a)s = {l~p)/2 



i-p 1 



^s{x, y) = x^y^ [Ci + C2 Inx] C^e''^ + C^e-''^ 



h)s^{l-p)/2 



tps{x,y) = x 2 y2 



Ci + C2X±V(l-P)'-4s2 



C^e''^ + C^e 



isy 



(3.8) 



(3.9) 



The solutions ( p.6| )-( p?9D are characterized by p and s as free parameters. Also, as was 
defined in ( p.5| ), variable y depends on arbitrary function a;(0). Bearing this in mind we can 
construct general solutions to the WDW equation by superposition of the given solutions. 
Later on we obtain the wormholes solutions for /c = 1. 



B. Case A{y) = Aq 



Now, we study the WDW equation with A{y) = Aq, where Aq is a constant, this case is 
equivalent to say that the effective cosmological term A(0) is proportional to the scalar field 
0, thus from Eq. (^.4|) we have the following wave equation 



(9^ f) 

X 7— r + pX- kx + AqX 

OX"^ OX 






i'{x,y) = 0. 



Setting ip{x,y) = e*^^X(a;), the equation for the variable x is 

x'^X"{x) + pxX' + (s^ + Aqx^ - kx^)X = 0, 
once again we will consider different values of k in this equation. 
t) k = ±l 



(3.10) 



(3.11) 



For the closed and hyperbolic model, the above equation is solvable by power series ||12[ , but 
the solution is a more complicated, then in order to look for another expressed in terms of 
standard functions 



we take X(x) = z f(z) 



X ), thus Eq. ( 3.11 ) transform into 



4z^f"{z) + {2+p + 8A)f + [aA^ + {p- 2) A + s^ ^ z^ + A^z^] / = 0, (3.12) 

in order to get an exact solution of this equation, we need to satisfy 2 + p + 8A = and 



4^2 + (p - 2)A + s^ = this imply that A = -^ and p = ~^^^f^ then, Eq. ( PIT^ 
reduce to a more simple form 



f"{z) + -{A,zTl)f = 0, 

next, we use the substitution f{z) = U{ip) with (2Aoipy 
wavefunction in terms of Airy functions is 



(3.13) 
AqZ =f 1, thus the universe 



V'(x, y) = {CiAt [- (2Ao)"^ (AqX^ t l) 



CoBi 



(2Ao)-t(Aoa;2Tl)]} 



X 

ii) k = Q 

In a flat FRW universe, the wavefunction is 

/A 



isy\ 



(3.14) 



ips{.x,y) =x -^"y 2" 



Ci/ 



M 



\^ 



+ C2I- 



M 



\^ 



c^e'^ + C^e 



■isy 



(3.15) 



where 1±m is a modified Bessel function and M = l\ i^-^) — s^, and by superposition we 
can obtain wavefunctions satisfying particular boundary conditions. On the other hand, for 
the case A; = ±1 we do not have a complete set of solutions, but only two particular solutions 
for the specific choice of the factor ordering parameter p. 



IV. GAUGE N=l/X 



In this section our aim is show another solvable quantum model with N = 1/x, then the 
action (Q.6\) becomes 



S 



1 



X x' — X y' + X A{y) — k 



dr. 



and with the following change of independent variables 

a = x'^ cosh(2?/), P = x'^ smh.{2y), 
the action ( [4 .11 ) takes the symmetric form 



S=\j \ ("" - /5") + Aia + A2/3 - fc 
where we have chosen the cosmological term as 



dr. 



A{y) = Ai cosh(2y) + A2 sinh(2|/), 
where Ai and A2 are constants, and the canonical momenta are given by 

P' 



TTn 



a 



^(3 



The Hamiltonian is constructed as 



1 



n = 2ni- 2n' + -{k - Aia - Aa/?) 



'^ ■ 2 
thus the corresponding WDW equation is 

92 



aa2 Qf]2 4 



52 1 

+ - (Aia + A2/3 - k) 



iP{a,P) = 0. 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



In particular, we can get exact solutions setting ip{a,P) = A(a)B(P) and introducing the 
new variables 



7 = (2Ai)-3(A; - As^ - Aia), r] = (2Ai)-3(A2/3 - 4^^), 



(4.8) 



where s is the separation constant. Then we obtain Airy equations G"{'y) — 76* = and 
H"{ri) — rjH = 0. Consequently, the universe wavefunction for this model is 



^P{a,(3) = \CiAi (2Ai)-3(A;-4s2-Aia) 



CoBi 



(2Ai)-3(/c-4s2-Aia) 



C.Ai 



(2A2)-3(A2/3-4s') +C4Bi (2A2)-3(A2/3 - 4s 



(4.9) 



To have a wavefunction for a particular theory we take the choice uj{(f)) = uq = const., we 
obtain a cosmological term in the Brans-Dicke theory that could lead to inflation in the 



classical theory |14,15 



A(0) = Ai0^+^ + \2(p- 



P+i 



where p is given by p 



a 
"=2 



{2uJo + 3)/3 then, variables a and (3 becomes to 



+ 



^-i 



in this case of Brans-Dicke theory, the universe wavefunction is 

^(a, 0) = iCiAi [(8Ai)-i {2k - 8s^ - A^a (^''+5 + 0"''+^)) 
+C2Bi [(8Ai)"i (2k - Ss^ - A^a (. 



hP+9 



C^Ai 
CiBi 



(8A2)-i (Asa (0"+^ + r^+^) - 8s') 
(8A2)-i (Aaa (0^+^ + 0- 



■p+tl_8s2 



(4.10) 



(4.11) 



(4.12) 



Depending on the sign of the argument the Airy functions are oscillatory or exponential, by 
superposition, different boundary conditions can be satisfled. In the next section we show 
some classical solutions for the general and some particular theories. 



V. CLASSICAL SOLUTIONS 

The classical equations of motion for this model derived from the action ([4.3|) are 

a"(r) - 2Ai = 0, (5.1) 



/3"(r) + 2A2 = 0, 
with the Hamiltonian constraint equation given by 

a'2 - /5'2 - 4Aia - AA2I3 + 4A; = 0. 
The solution of system (|5rT|)-(|57^) is 

a(r) = AiT^ + CiT + Ca, 



/?(r) 



-AsT^ + C^T + C4, 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



where Ci and C2 are integration constants and satisfy Cf — C^ — 4A1C2 — 4A2C4 + Ak = 0. 
In the following we study simply choices of the uj{(j)) function, this lead us to Brans-Dicke 
and Barker theories, also we consider three general parameterized theories that have been 
studied by other authors 0,0 in the Appendix. 
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A. Brans-Dicke theory 



Choosing u;(0) = ujq, then y = | ln0 and we obtain the following solution 

a{T) = [(Ai + A^y + (Ci - C,)r + C^ - C4]™ 
"(Ai - A^y + (Ci + C^)t + C2 + C, """ "" 



(5.6) 



(Ai - A^y + (Ci + C^)t + C2 + C^ 
(Ai + A2)r2 + (Ci - C3)r + C2 - C4 



(5.7) 



where m = (1 + p)/4:p, n = l/2p and p = J{2ljq + 3)/3. Here we have a family of solutions, 
depending on several parameters, in the next subsection we take particular values of the 
parameter to have explicit solutions that are easier to analyze. We found singular and non 
singular solutions. 



1. Non singular solution 



Setting Ci = C3 = C4 = 0, C2 = k/Ai and Ai = A2 in the above equations, we have the 
solution 



a{T) = aQ 



+ 1 



and 



(r) 



i + 1 



(5.^ 



(5.9) 



where qq = (k'^/Af)* and Tq = k/2Al, this solution corresponds to a cosmological term of 
the form A{a,(3) = AiJ{a + {3) /{a — j3) and is not singular for /c = 1, we can verify this, 
by direct substitution of the solution in the metric invariants 



i?= -6 



+ 20^a'^ + a(t)a'(t)' + a(t)'^a" 



(5.10) 



^^-4^ 



3j,„/j,/ „3i2„// 



k — a (ha (h— ad) a 



D _ ^ d3/2 



(5.11) 
(5.12) 



i?3 



Rl 



(5.13) 



Introducing the solutions (|5.8|) and ( |5.9|) in the Ricci scalar ( |5.1CI|) we obtain 

6 



R 



2^4 / 1-2 I j^ 



anT, 



o'o 



2(m+l) 



Hr^ + n + 



+ 2t'{ kr^ + a^m(6m - 2n - 1) ^ + 1 )+2almT^ ^ + 1 



Tn 



.^0 



and the -Ri invariant takes the form 

3 



R^ 



Aatrl (^ + l) 



4(m+l) 



M< + r') + 



r 



r 



n2 



+ 2r' <^ A;ro' + 2a^m(2n - 2m + 1) ^ + 1 U - 2a^mro' ^ + 1 



It is easy to check that this two invariants are finite for all values of r. 



(5.14) 



(5.15) 



2. Singular solution 

Another classical cosmological solution can be obtained by choosing Ci = C3 = 0, C2 

C4 = k/2Ai and Ai = A2, thus the solutions ( ^.6|) -( p?7D are reduced to 



a{T) = ay"", 



(5.16) 



r 



/ -2n 



(5.17) 



1 — — 
where Oq = 2iTQ ^ and 0o = ''"o"- These solutions in the cosmic time t take the form 

2(l+p) 

a{t) = a't^^^, (5.18) 



0(t) = 0't3^ 



(5.19) 



2(1+P) 



where a' = a'^ {2^/^t^/^p (Ig^)} ^-^^ and 0' = 0[) {21/^3/4^ (^l^)}^''-^ For other theories 
see the Appendix. 



VI. QUANTUM WORMHOLE SOLUTIONS 



In this part we show some special solutions of the WDW equation for A^ 
A^ = 1/x, which are known as quantum wormholes. 



1 and 
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A. Case A(y) = 0, A^ = 1 
For the closed universe {k = 1), assuming p = 1 the solution to WDW equation (|3.5|) is 



'>Ps{x,y) 



CJ,s (^) + C,K,, (^' 



c^e'^y + ^46-^^2' 



(6.1) 



If we choose the family of solutions 



ilj,{x,y) = Ce''yK,, 



X 



(6.2) 



then we can generates a quantum wormhole basis, namely, making wave superposition 



tpwH{x,y) 



~oo / -y^ \ 

^ C{s)e''yK,J-\ ds, 



(6.3) 



if we take C{s) 
||T6ll we obtain 



^ifis ^^ _ Qonst.), by means of a Kontorovich-Lebedev integral transform 



^/.(a;,^) 



^-i- cosh [22/+^] 



(6.4) 



This solution is a quantum wormhole, because it satisfies the Hawking-Page regularity 
boundary conditions, i.e., the wavefunction is exponentially damped for large spatial geo- 
metry, and also, is regular when the spatial geometry degenerates [[7|. Wormholes may play 
an important role in solving problems associated with the complete evaporation of black 
holes, and it is believed that they produce effective interactions in the low energy physics 



17| that turn the coupling constants of nature into dynamical variables [|T8[. The set of 



wavefunctions ips and ip^ are spanning the same space of physical states, and are both bases 
of the Hilbert space of the model in the corresponding representation, and Eq. (|6.3| ) is the 
connection between these bases tps and ^^ [jl9|,^. For Euclidean wormholes in BW theory 



see 21 . 



B. Case k{y) = Kq, N = I 

In a similar way, for a flat model {k = 0), we find for Aq < the following quantum 
wormhole basis 



'^t,{x,y) 



/Aqx^ cosh [3y—^] 



(6.5) 



again this solution satisfy the mentioned Hawking-Page wormhole condition, and by direct 
substitution we can verify that basis (|6.5|) is a particular solution of the Eq. (p.lO| ). 
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C. Case A{rj/0 = 0, N = l/x 

Also quantum wormholes can be obtained from action ( |2.6| ) with N = l/x, if we introduce 
the following independent variables 



^ = X cosh y, 



7] = X sinh y, 



then the new action is 
1 



*=2 



ee - ev" - vv + vv - «= + (f - vrM.„/o 



dr. 



and the corresponding WDW equation is 



^.-^.-(^'-v')k+{e-vrAiv/o 



H^,v) = 0- 



(6.6) 



(6.7) 



(6.8) 



Choosing A{ri/C,) = and considering a closed model {k = 1), Eq. ( |6.8| ) simplifies to 



532 Q2 

— - e' - — + ^' 

^^2 ^^2 



v^(e,^) = o, 



(6.9) 



we can see that this is, the equation for two harmonic oscillators with opposite signs of the 
energy with solution 



'^{^,v)wH = Hn{i)Hn{'n)e 



(6.10) 



where H^ are Hermite polynomials of order n. These harmonic-oscillator solutions are a 
quantum wormhole basis, since satisfied the Hawking-Page regularity boundary condition, 
(i) the solutions are regular at the origin and (ii) damped at infinity. 

VII. THIRD QUANTIZATION 

It is known that the WDW equation is a result of the quantization of a geometry and 
mater (second quantization of gravity). The procedure of quantizing the wavefunction ?/> 



is called third quantization [11], in this theory we consider ip as an operator acting on the 
state vectors of a system of universes and can be decomposed as 



?/'(«,/?) 



ds 



C(s)^,(a,/3) + C(s)V*(«,/5) 



(7.1) 



where "0(0, (3) and ip*{a, P) form complete orthonormal sets of solutions to the WDW equa- 
tion. This is in analogy with the quantum field theory, where C{s) and C^{s) are creation 
and annihilation operators. Thus, we expect that the vacuum state in a third quantized 
theory is unstable and creation of universes from the initial vacuum state takes place. In 
this view, the variable a plays the role of time, and variable /3 the role of space. The 
wavefunction ^(a,/?) is interpreted as a quantum field in minisuperspace. 
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We assume that the creation and annihilation operators of universes obey the standard 
commutation relations 



[C{s),C\s')]=6{ 



s — s 



[C{s),C{s')] = [C\s),C^{s')] = 0, 



(7.2) 



and we expand the field ip in normal modes ips as shown in Eq. ( [7.11 ). Here, the wave number 
s is the momentum in Planck units and is very small. Let us consider now the last quantum 
model ( ]4.7| ) for A2 = and closed universe k = 1. Then, the WDW equation becomes 

-(Aia- 1) 



da^ 9/32 4 



i/j{a,f3) =0. 



(7.3) 



A complete set of normalized positive frequency solutions to this equation ||15 

is/3 



IS 



^+(a,/3) 



(I6A1 



Ai 



iBi 



(2Ai)"3(i-4s2- Aia) 



(2Ai)~3(l-4s^- Aia) 



+ 



(7.4) 



and 






(16A 



(2Ai)-3(i-4s2-Aia) 



+ -e 3Ai Qi 
2 



(2Ai 



;i-4s2-Aia) 



(7.5) 



the wavefunctions '0^ and ip~ can be seen as a positive frequency out going and in modes, 
respectively, and these solutions are orthonormal with respect to the Klein-Gordon scalar 
product 

(V^., V^s') =^J^s^(, rs'd/3 = 5(3 - s'). (7.6) 

As both sets, ( [7. 41 ) and ( [7.5] ) are complete, they are related to each other by the Bogoliubov 
transformation defined by 



^+(a,/5) 



dr 



M+{s, r)^-(a, p) + M-{s, r)V',"(a, /3) 



(7.7) 



The wavefunction ilj{a,f3) can be expanded in terms of annihilation and creation operators 
for each set of modes, then with wavefunctions ( [7.4| ) and ( [7.5| ) we obtain that the Bogoliubov 
coefficients M~^{s, r) = 6{s — r)v{s) and M~{s, r) = S{s + r)w{s) are 



M^{s,r) = 6{sTr)-^ 



1 (1-4^2)^ 

2 



(l-4s^)^ 

+ e 3Ai 



(7.8) 



The coefficients M^{s,r) are not equal to zero. Thus, two Fock spaces constructed with 
the help of the modes ( [7.4| ) and ( |7.5| ) are not equivalent and we have two different third 
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quantized vacuum states (voids): The in- vacuum | 0,in) and out- vacuum | O,out), (which 
are the states with no FRW-hke universes) defined by 



Cin{s) I 0,in) = and Cout{s) \ O,out) = 0, 



(7.9) 



where s G R. As was mentioned, since the vacuum states | 0,in) and | O,out) are not 
equivalent, the birth of the universes from nothing may have place, where nothing is the 
vacuum state | 0, in). The average number of universes produced from nothing, in the s-th 
mode N{s) is 

(7.10) 



N{s) = (0,m I ClUs)Coutis) I 0,m 
as follows from Eq. ( [7. 81) we get 



N{s) 



1 (l-4s^)? (1-43^)^ 

_g 3Ai g 3Ai 



(7.11) 



In view of Coleman's wormhole mechanism [18| for the vanishing cosmological constant, we 



take the constraint Ai < l/Svr x 10 ^^^m^, also assuming \s\ << 1, then notice that we can 
get the number of state N{s) 



Nis] 



1 2 (l_4s2)i 

r-^ p 3A]^ ^ '' 



(7.12) 



This result from third quantization shows that the number of the universes produced from 
nothing is exponentially large. 

VIII. UNCERTAINTY RELATION 



In this section we will study third quantization for the gauge A^ = 1, in order to get the 
uncertainty Heisenberg relation, we start from WDW equation ( |3.4| ), setting A{y) = Aq we 
have 



X 



dx' 



+ px 



d d^ 



dx dy"^ 



X\k - AqX^) 



^(x,y) = 0, 



(8.1) 



if we choose p = —1 and introducing the new variable z = x'^, the above equation transform 
into 



d^ 



1 92 



dz"^ 4^2 dy'^ 



+ Uiz) 



'^{z,y) = 0, 



^.2) 



U{z) = --{k-Aoz). 



(8.3) 
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The third-quantized action to construct de WDW equation 



IS 



-^ 



l)^-i(|)^--w^'^' 



dzdy. 



•4) 



Now we use a Fourier decomposition in order to obtain decoupled degrees of freedom, the 
universe field ip{z, y) will be expand in sine and cosine functions. We assume that our 
quantum system is confined to a one-dimensional box with periodic boundary conditions, 
with the coordinate length fixed at a arbitrary value M, then 



.5) 



^iz,y) = —— lij{z,0)+ Yl -^bP+iz,q)cosqy + ^-{z,q)smqy] 

^^^ y q=2^n/M V 2 

where ^J+i^z, —q) = ip^{z^q) and ^J-i^z^ —q) = ip_{z,q). By substituting the expansion (^]5| 
into Eq. (^.41) , we obtain that 






'^)^-^*'--(^« 



(8.6) 



where we have denoted the mode variables ^(-2,0) and "^±{z^q) by V'(t(^) and Eq. (| 
is rescaled to '^^ — 



M/2'Kil)cj. The above sum includes zero mode ip{z, 0) for each pair 
(g, — g), in this way the mode variable ■j/'o- is completely decoupled from each other. Now the 
Hamiltonian is given by 



h.=y: 



^l+[\ + U{z)\^^l 



(8.7) 



To quantize this system (|8.7| ), we impose the equal time commutation relations 

i)a,K =iSa,a'- (8.8) 

The wave functional of any energy eigenstate is factorized as 

v|/ = n^4^,^a]- (8.9) 



In order to get the Heisenberg uncertainty relation we use the Schrodinger picture, then we 
make the substitution ip —^ ip^ and ir^ -^ ~^'A 
mode variable is 



— z-7^. Then, the Schrodinger equation for each 



Ma 

dz 



1 



d^ 



di)l 



4^2 



+ U{z) ^l 



^. 



(8.10) 



We will solve the above wave functional equation by using the Gaussian ansatz 
'ilJa[z,^pa] = C exp{--Aa{z, q)['ilj^{z, q) -e^{z,q)f + iP^{z,q)['ilj„{z,q) -e^{z,q)]}, (8.11) 
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Aaiz, q) = D„{z, q) + iF„{z, q), 



(8.12) 



where the real functions D„{z, q), F„{z, q) Pa{z, q), and 6*0- (z, q) have to be determined from 
Eq. (^.101) . C is a normahzation of the wave function. The inner product of two functional 

ipiaizjipo-) and 4'2a{z,ip(T) is defined by 



< ^1 I ^2 >z= / dllJaiJla[z,tp^]tp2a[^,tp^]. 



(8.13) 



We shall calculate Heisenberg's uncertainty relation, the dispersion of ipa is given by 
{Ail)„f =< 4jI>,- < ip^ >2, from Eqs. ( pTD and (|1|) we have 



<€>. 



2D^{z,q) 



+ Vliz,q), <i)a>z=vliz,q), 



then 



(A^.) 



2D^{z,qy 



(8.14) 



(8.15) 



and the dispersion of n^ is (AvTo-) =< vr^ >2 — < n„ >l, where 



< < >-. 



D^{z,q) , F^{z,q) 



2D,{z,q) 



+ P^{z,q), <n„>,= P^{z,q), 



(8.16) 



then simplifying 



(Att^) = h 



thus we obtain the uncertainty relation 



1 



{Ai;^)'iAn^)' = -il 



2D^{z,q) 



Fa{z,q) 
Da{z,q) 



(8.17) 



In order to obtain F„{z, q) and D„{z, q) we substitute the ansatz ( B.lll ) into Eq. (|8 
we obtain the equation of motion for ^^(z, q): 

d q^ 

- ^j-^K{z, q) = -Al{z, q) + j^ + U{z). 



(8.18) 
I|), thus 



Now we substitute 



A„{z,q) = -i—lnu^{z,q), 



^.19) 



(8.20) 



where Uo-{z, q) is the solution of the WDW equation 

Ua{z,q) = 0. 



92 q2 



Q^2 42;2 



(8.211 
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Now first, we analyze the behavior of u^ for large scales (a -^ oo), then the second term 
q^ /Az^ in the wave equation ( p.21| ) can be neglected, and its solution in terms of Airy 
functions is 



Ua{z,q) = Ai 



(aAI) ' {Aoz - k)] + gBi UaAI) ' (Aqz - k) 



where ^ is a complex constant. Substituting (|8.22| ) into (|8.20|) , we obtain that 



2\i 



Ai 


(4A2)-^Ao^-A;) 


+ gBi' 


(4A2)-^ (Ao^ - A;) 


Ai 


(4A2)-i (Ao^ - k) 


+ gBi 


(4A2)-^ (Ao^ - k) 



A,(z,g) = zAo(4A^)3 



where the prime denotes differentiation with respect to z, also we have 

L(Im^)(Ai' + ^Bi') 



D,{z, q) = ReA^{z, q) = Ao(4A; 



and 



F^{z,q) = -lmA^{z,q) = Ao(4A, 



therefore the uncertainty relation becomes 



I Ai + ^Bi |2 ' 
2.1 Ai'Ai+ I ^ |2 Bi'Bi + (Re^)(Bi'Ai + Ai'Bi) 



Ai + ^Bi 



(A^^f (Att.) 



1 



n 



Ai'Ai+ I ^ p Bi'Bi + (Re^)(Bi'Ai + Ai'Bi) 



4 [ (Im^)2 
For large scales (a — > oo), the above uncertainty relation is given by 



(AV.)' (Att,,) 



1\ 1 

4\ (lm^)2 



(IH'-l)sin 



^{4AlrHk-Aoz)l + ^ 



2(Re^) cos 



liAAl)-Hk-Aoz)l + l 



for g = —i, we get that 



{A^^f (Att^) 



and when g ^ —i and | /? |>> 1 



(A^.)'(A7r,) 



4 



1 + 0(1/3 



(8.22) 



(8.23) 



(8.24) 



(8.25) 



(8.26) 



(8.27) 



(8.28) 



(8.29) 



This means that the quantum fluctuations of the third-quantized closed universe field are 
bounded at a finite value according Eq. ( ^.26| ), during the universe expansion (a -^ oo). 
Now we study the behavior of u^ for small scales (a — ;> 0). Assuming that < Aq < 1 in 
Planck units and k = 1, the solution to Eq. ( p. 211) is 
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Ua{z,q) 



/.(|)+7A-.g 



(8.30) 



where Ii, and Ki, are modified Bessel functions, 7 is a complex constant, and u = \V^ " Q^- 
Substituting tlie general solution ( |8.3CI| ) into Eq. ( ^.201 ), we obtain 



D. 



ReAa- 






(8.31) 



ImAn 



Z-' I h + iK P +I'Ju+ I 7 P KKu + Rei{Kh + KK,) 



(8.32) 



(8.33) 



2 I /, + 7K, |2 
tlie prime denotes differentiation witli respect to z. Tlien, tlie uncertainty relation is 

where 

a = {Im^Y^z^ [z-\ll+ I 7 I' Kl) + /:/,+ I 7 p <K, + Re7«/. + ^K,)] , (8.34) 

since we have assumed that our system is confined to a one- dimensional box with periodic 
boundary conditions, with the coordinate length fixed at an arbitrary value M, if we take 
the limit M -^ 00, then we can choose z/^ > 0, and for small scales (a — > 0) the asymptotic 
behavior ||13[ of Eq. ([OBI ) is 



(A^<,)2(A7r.)2 = l{l + gz-4H}^ 



(8.35) 



where Q is some positive constant. This means that the fluctuation of the third-quantized 
universe field becomes large for small scales (a ^ 0). For the case in which z/^ < 0, so that 
V is pure imaginary, we obtained an oscillatory behavior of G^ and does not have a definite 



magnitude |22| when a ^ 0. 



IX. CONCLUSIONS AND DISCUSSION 

In this paper we have investigated some issues of quantum cosmology of simple mini- 
superspace models in the Bergmann- Wagoner theory, after quantizing the models we have 
obtained a class of separable Wheeler-DeWitt equations. We construct universe wavefunc- 
tions for fiat, closed and open cosmologies. 

The first physical consideration of wavefunctions here obtained, is related with the fact 
of that in quantum gravity the topology of spacetime is expected to fluctuate on Planck 



scales I^Bj. The spacetime might be multiply connected, these connections are microscopic 
wormholes. Hawking and Page ^ have argued that wormholes are solutions of the WDW 
equation, which are exponentially damped for large three-geometries, and regular in some 



1^ 



suitable way when the three-geometry collapses to zero. We have shown two discrete basis 
of wormholes for A(0) = and A(0) = 3Ao0 cosmological terms. Also, we have presented 
singular and non singular general solutions in the closed Brans-Dicke cosmology. 

The second physical consideration is related with the third quantization of the universe. 
The solutions of WDW equation have problems with the probabilistic interpretation. In the 
usual formulation of quantum mechanics a conserved positive-definite probability density is 
required for a consistent interpretation of the physical properties of a given system, and the 
universe in the quantum cosmology perspective, does not satisfy this requirement, because 
the WDW equation is a hyperbolic second order differential equation, there is no conserved 
positive-definite probability density as in the case of the Klein-Gordon equation, an alter- 
native to this, is to regard the wavefunction as a quantum field in minisuperspace rather 
than a state amplitude, and the strategy is to perform a third quantization in analogy with 
the quantum field theory. Thus, the aim of this approach is to construct a consistent pro- 
babilistic measure in quantum gravity by promoting the wavefunction of the universe to a 
quantum field operator that acts on a Hilbert space of states. The vacuum state in this space 
is identified as the state where the universe does not exist. Topology changing processes 
can then be described by including self-interaction of the universe field. In quantum field 
theory, the particles are created from the vacuum by a time-varying external potential and 
this suggests that universes could be created via a similar process. In the third quantization 
the universe field is expanded into positive frequency in- and out-mode functions and their 
hermitian conjugates. The in- and out-modes are related to one another by the Bogoliubov 
coefficients and these determine the number of universes in a given mode. We have found for 
a simple minisuperspace model with a potential term X{y) = Ai cosh {2y), that the number 
of the universes produced from nothing is exponentially large. 

On the other hand, quantum fluctuations of flat and open minisuperspace models have 
been analyzed by the method of a time dependent Hermitian invariant [^, and was ob- 
tained that the fluctuation around each WDW trajectory converges to its minimum rapidly 
in the course of the cosmic expansion. Nevertheless, the closed models have not been ana- 
lyzed by this method. We studied the third-quantization of a closed minisuperspace model 
by using the functional Schrodinger equation in order to investigate the Heisenberg uncer- 
tainty relation, and we found that quantum fluctuation of the third-quantized universe fleld 
becomes small for large scales (a — »• oo) in the course of cosmic expansion. Also we found 
an exponentially large dominance of quantum effects for small scales (a — ^ 0). The opposite 



behavior for the uncertainties was obtained by Phole |25], the difference is do to the fact 



that he was considering the classically forbidden region as was pointed out by Horiguchi 



22| 
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APPENDIX A: CLASSICAL SOLUTIONS FOR SOME SCALAR-TENSOR 

THEORIES 



Here we consider some particular choices for the couphng function uj{(f)), that have been 
used previously by some other authors. A common property of these coupling functions is 
that they give an analytical result for the integral in Eq. (|2.5| ). 



1. Barker's theory 

The choice ^(0) = |^^ lead us to the Barker's theory. This coupling function has the 



20-2 

characteristic of making G = to first order in the weak-field limit 
function we obtain the solution 



and 



r(Ai + A2)r2 


+ (Ci- 


-Cs)r 


+ C2 


-c. 


[(Ai-A2)r2 


+ {Ci + Cs)t 


+ C2 


+ C4 



a{T) = < tan In 

i LI- 

[[(Ai + A2y + (Ci - Cs)t + C2- C4] 
[(Ai - A^y + (Ci + C3)r + C2 + C4]]i 



V3 



(t) = tan In 



r(Ai+A2)r2 


+ (Ci- 


-C,)t 


+ C2 


-c. 


_(Ai-A2)r2 


+ (Ci + C,)t 


+ C2 


+ C4 



V3 



Using this coupling 



+ 1. 



(AI) 



(A2) 



2. Theory 1 

For the coupling function 2ci;(</)) + 3 = Bl(j?'^^^^\ i?i > [^], we obtain the following 
solution 



a[T) 



In 



•(Ai + A2V + {Ci - C^)t + C2 - C4 
(Ai - A2)r2 + (Ci + C^)t + C2 + C4 

[[(Ai + A2)r2 + (Ci - C^)r + C2- C,] 
[(Ai - A2)r2 + (Ci + C,)r + C2 + CJ]^ 



V3(i + 1) "\ 2(1 + 1) 



(A3) 



r 



In 



(Ai + A2)r^ + (Ci - C^)t + C2- C4 
(Ai - A2)r2 + (Ci + C3)r + C2 + C4 



V3(i+1) "> ! + l 



(A4) 
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If we take 2a;(0) + 3 
equations is given by 



3. Theory 2 

B2 I In ( ^j 1^^"^ IJlOl, the resulting solution to the cosmological 



a(^) = 00 ' exp 



Hvi'^-'"" 



a. (5/1 










[(Ai + AsK^ 


+ (Ci- 


-C3)r 


+ C2 


-C4 


[(Ai-A2)r2 


+ (Ci + C3)r 


+ C2 


+ C4 



1 

l-<5 



{[(Ai + A^y + (Ci - Cs)t + C2- C4][(Ai - A^y + (Ci + C3)r + C^ + C,]}' , (A5) 



and 




(r) =0oexp<^ W^(l-(5) In 



(Ai + A2)r2 + (Ci - C3)r + C2 - C4 



(Ai - A2)r2 + (Ci + C^)t + C2 + C^ 



1 
1-4 



(A6) 



6 = 1 



a[T) = <pQ exp 



(Ai + A2)r^ + (Ci - C^)r + C2 - Q 
(Ai - A2)r2 + (Ci + C3)r + C2 + Q 



{[(Ai + A2)r2 + (Ci - C3)r + C2 - C4][(Ai - A2)r2 + (Ci + C^)r + C2 + Q]}' , (AT) 



and 



[.T) = 00 exp 



(Ai + A2)r2 + (Ci - Cs)t + C2- C4 
(Ai - A2)r2 + (Ci + C3)t + C2 + C4 



1 



4V3B2 



(A8) 



4. Theory 3 

For the last function considered here we choose 2c(j(0) + 3 = — ^'^ . ^ with A > 0, B^ > 
lOH, and we obtain 



1-1 It 



a(r) = (pf)'^ ll — tanh In 



(Ai + A2)r2 + (Ci - C3)r + C2 - C4 
(Ai - A2)r2 + (Ci + C3)r + C2 + C4 



-, A JL 



1 

" 2A 



{[(Ai + A2)r2 + (Ci - C,)r + C^ - CM^i " ^2y + (Ci + C3)r + ^2 + ^4]}% (A9) 



and 



(r) = 00 S 1 — tanh In 



(Ai + A2)r2 + {C, - C3)r + C2 - C^ 
(Ai - A2)r2 + (Ci + C3)r + C2 + C4 



A /3 1 A 



(AlO) 
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